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I. Introduction and main results. 

There exits an increasing literature devoted to the study of the killed Brownian motion 
on unbounded domains and the associated survival probability, for instance see [C.Z.], 
[K.], [L.], [P.]. Here, we consider the large time asymptotic of the heat kernel in general 
Benedicks domains with Dirichlet boundary conditions. The case of Benedicks domains 
with compact holes was treated in [C.M.SM.l]. Although some cases with non compact 
holes can be studied using the same methods, it seems hopeless to extend this approach 
to the general case. In the present paper we derive a new inequality (of Nash type) which 
allows to control the general case of Benedicks domains. For these domains we answer 
positively a conjecture of Davies [Dl.] about limits of ratios. The case where there is 
only one positive harmonic function vanishing on the boundary (up to a positive factor) 
was treated before for more general domains in [C.M.SM.l] and for exterior domains in 
[C.M.SM.2]. See also [B.J.]. The most interesting case is of course the case when the cone of 
non negative harmonic functions with Dirichlet boundary condition is of dimension larger 
than one. We will see that for Benedicks domains, the asymptotic rate of convergence to 
zero of the heat kernel (or of the survival probability) clearly reflects this property. 

We recall that a Benedicks domain O is the complement of a closed set D c x {0} where 
D is an open non trivial subset of the hyperplane xa = (see [B.]). We also assume as in 
[B.] that each point of D c x {0} is regular for the Dirichlet problem in O. If x e R d we 
will sometimes write x = (x, Xd) with x G R . We recall that for Benedicks domains the 
cone Vq of non negative harmonic functions with zero boundary condition is of dimension 
one or two (see [B.]). An integral test for this dimension is given in [B.] (see also [A.]) while 
Theorem 1 below gives a probabilistic characterization. In the case of dimension two, we 
denote by u\ the extremal of the cone normalized by 

lim 

We recall that since u\ is an extremal we have 

lim = . 

x d ^-oo Xd 

We denote by u 2 the symmetric function with respect to the hyperplane xa = 0, namely 

u 2 {x,x d ) = iti(x, -x d ). 
It is another extremal function which satisfies 

lim ^A = Q and l im =u 

x d ^+oo Xd x d ^-oo Xd 

Note that u\ and u 2 generate Vq. We also denote by v s the function 

v s (x, x d ) = ui (x, x d ) + u 2 (x, x d ) , 
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which is modulo a constant factor the unique symmetric element in Vq. 

We denote by pt(x,y) the heat kernel of Q with Dirichlet boundary conditions. For 
a Brownian motion X t in O with absorption at the boundary, let Tq denote the exit time 
from O. We will denote this random variable by T when there is no ambiguity on the 
domain. P X (T > t) will denote the probability of survival up to time t starting at the 
point jGD, namely 

P x (T>t) = [ p t (x,y)dy. 
Jo, 

We will also use the heat kernel of the half space xa > denoted by pf(x,y). We now 
formulate our result. 

Theorem 1. For any Benedicks domain O, for each (x, y) G O x O, the function of time 
t 1+d / 2 p t (x 7 y) has a limit when t tends to infinity. 

If this limit is infinite for one pair (x, y), it is infinite for any pair (x' , y') G O x O, the cone 
Vq, is of dimension one, and we have for any x±, • • • , x^ in Q 

li m Pt( x i> x z) _ u(x 1 )u(x 2 ) 
t^oo Pt (x 3 , X4) u{xz)u{x4) 
where u is any function in the cone. 

If the limit is finite for one pair (x,y) (and hence for all pairs of points), the cone Vq has 
dimension two and we have for any x and y in Q 

2 

l}n^t 1+d/2 pt(x,y) = ( u i( x ) u i(y) + u 2 {x)u 2 {y)) . 

In the first case, the exact asymptotic rate in time of pt{x, y) differs according to the 
domain. For example, in dimension two, for the exterior of a segment one gets t -1 (logt) -2 , 
and for the complement of a half line one gets t~ 3 / 2 (see [C.M.SM.2] and [B.S.]). 

A result similar to Theorem 1 holds for the survival probability. 

Theorem 2. For any Benedicks domain O, for each x G O, the function of time t^P^ (T > 
t) has a limit when t tends to infinity. 

If this limit is infinite for one point iGfi, it is infinite for any point j/GO. The cone Vq 
is of dimension one, and we have for any x, and y in O 

.. P x (T>t) u(x) 

hm ; = , - 

t^P y (T>t) u{y) 

where u is any function in the cone. 

If the limit is finite (for one point), the cone Vq, has dimension two and we have for any 
x G O 



lim t 1 / 2 P x (T>t) = 



v s (x) 



Here also, in the first case the asymptotic rate depends on the domain. 

The rest of the paper is organized as follows. In section II we prove Theorem 1. The 
proof of Theorem 2 is given in section III. In the appendix we give a proof of symmetry 
estimates needed in the proof of Lemma 4. 
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II. Asymptotic of the heat kernel. 



The proof of Theorem 1 will be a consequence of several results. We first establish a general 
estimate which may be of independent interest. This Lemma follows by applying Nash 
inequality to the heat kernel. In fact, an analytic proof is based on the inequality \\4>\\\ < 
C2IIV0IIIH0II?, for 4>{x) = p t (x,y). A direct computation yields \\p t (-,y)\\i = P2t{v,v), 
l|V.Pt(-,2/)||! = ~P2t(y,y) and ||pt(-,y)||i = P y (T > t), leading to an inequality of the type 
(1) below, with p^ t instead of pzt- In the sequel we give a direct and simpler probabilistic 
proof. 

Lemma 3. For any x and y in and any t > 0, we have 

P3t(x,y) < — -^P X (T > t)Py(T > t) . (1) 



Proof. By the semi group property we have 

P3t(x,y)= / dz\ / dz 2 Pt(x,z 1 )p t (z 1 ,z 2 )p t (z2,y) ■ 
Jo, Jn 

On the other hand from the Gaussian bound 

v .(x v) < - e -\ x ~y\ 2 / 2t 

we obtain 

Pt{zuz 2 ) < jA^e- ]zi - Z2]2/2t < 1 



(2yrt) d / 2 - (27vt) d / 2 ' 

The result follows from the above expression for P X (T > t). QED 

By a direct computation, using the results of [B.S.] or [D2.], one can check that for 
cones the following limit exists 

lim ( 2 * t ) d/2 Pt(xi,X2) 



t-oo P Xl (T>t)P X2 (T>t) 

This limit also exists in the case of exterior domains as follows from the results in [C.M.SM.2]. 
For the case of Benedicks domains with two extremal harmonic functions this follows from 
Theorems 1 and 2. 

Lemma 4. For any Benedicks domain Q with a cone Vq of dimension two, there is a 
constant F > 1 such that for any t > 1 and any (x, y) G O x O we have 

Pt(x, y) < (1 + v s (x))(l + v s (y)) (2) 

We also have for some finite constant K, for any t > and any x e 0, the bound 

pat > t) < kI±$L . 
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Proof. We first observe that if V > 4/(2-7r) d / 2 , the first estimate is true for 1 < t < 4 due 
to the Gaussian bound and the non negativity of v s . We are now going to prove recursively 
that if (2) holds on the time interval [l,3 n ] (where n is a positive integer) it also holds 
with the same constant on the time interval [3 n , 3 n+1 ]. 

We will assume that yd > 0, the case yd < is similar, and yd = is recovered by 
continuity. We first prove an upper bound for Py(T > t), using the recursive bound (2) 
(i.e. assuming t E [l,3 n ]). Let a = r" 1 /^). We have 

Py(T>t)=h+I 2 , 

where 

h= Pt(z,y)dz and I 2 = I p t (z,y)dz . 

J\z d \<aVt J\z d \>aVt 

We first estimate I 2 . Recall that (see [B.]) v s (z) > \zd\. Therefore 

h = f ]~^ Vs Z p t (z,y)dz < 1 / (l + v s (z))p t (z,y)dz 

J\z d \> a Vt 1 + v s(Z) 1 + aVt J\z d \>aVi 

< ~r~~~~ ~n /(I + Vs(z)) Pt (z, y)dz < ±±^M . 

1 + OL\Jt J (X\Jt 

In order to estimate I\ we will split this integral into two parts. Namely 

h = Ji + J 2 

where 

Ji = / p t (z,y)dz and J 2 = _ p t (z,y)dz . 

J\z d \<aVt , \z-y\<aVd-lVt J\z d \<aV~t , \z-y\>a^J d-lVt 

We first estimate J 2 - We observe that 

d-l 

J 2 <^J 2 ,z with J2,i= p t (z,y)dz . 

1=1 J\z d \<aVt , \zi-yi\>aVt 

We now estimate each J 2 i in terms of I 2 . To this end we use the symmetry with 
respect to the hyperplanes Zd = ±(^z — yi) denoted by . Using Lemma A from the 
appendix we get for zi — y\ < —\zd\ 

Pt(z, y) < p t (Sfz, y) + p t (Sfz, y) 

and a similar estimate if zi — y\ > \zd\- Therefore 

J 2 ,i < 4/ 2 
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and we get a bound by the previous estimate of I^- 
We now estimate J\. We have 



Ji = / _ p t (z,y)dz = 

J\z d \<aVt , \z-y\<c i Vd=TVt 
J\z d \< a Vt ,\z-y\< a Vd^TVt {1 + V S {Z)) 1 / Z 



and using Schwartz inequality we get 

Ji < 

/ _ p t {z,y)(l + v s (z))dz I 

\J\z d \<aVt , \z-y\<aVd=TVt J\z d \<aVt , \z-y\<a^/dTl 



<(l + v,(y)) 1/2 (f *k 

\J\z d \<aVt , \z-y\<aVd^TVt I 1 + V 



_Pt(z,y) 

x 1/2 

Pt(z,y) 



1/2 



1/2 



dz 



(z)) dZ 



In the last integral we use the recursive assumption (2) and get 

a < (i + vM)^t-w^ ( I 7t£t7&* 

\J\za\KaVt , \z-y\<aVd^lVt \ l +V S {Z)) ^ 

<C d (l + v s (y))r 1/2 a d/2 t- 1/2 , 

where Ca is a constant depending only on the dimension d. 

Collecting the estimates and using the relation a = T~ 1 ^ 2+d \ we get for a positive 
constant C4 (which can always be assumed larger than 1) independent of T, x, y and 

te[l,3 n ] 



(T > t) < C4(1 + ^ (y)) (1 + a^F 1 / 2 ] < 2 C 4(l + ^(y)) r i/(2+d) 



This is the second part of the Lemma for t G [1, 3 n ] if we take 

K = 2C 4 r 1 /(2+^) . 

For t G [0, 1] the second part of the Lemma follows at once if we take K > 1 since 
P y (T > t) < 1. Using Lemma 3 we get for t G [3 n , 3 n+1 ] 

If we now chose 

(2C 4 ) (4+2d)/d 



r > 



( 27r )(2+d)/2 
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we obtain for t G [3 n , 3 n+1 ] the required estimate, and we complete recursively the proof 
of the first part of Lemma 4. QED 

We now come to the proof of Theorem 1. First of all, we observe that if there is a pair 
(x,y) G fi x fi such that v- 1+d l 2 'p t {x,y) has a finite accumulation point when t tends to 
infinity, then using Harnack inequality we conclude that t^ l+d / 2S> p t {x, x) has also a finite 
accumulation point. This implies as in [C.M.SM.l] that we can find a diverging sequence 
(t n ) such that for any s > and any (y, z) G fi x fi, the sequence (t n + sp 1+d ^ 2 'pt n + s (y 1 z) 
converges to a function f(y, z) which is harmonic in y and z (separately) and vanishes at 
the boundary of fi. In the case of manifolds, see [A.B.J.] for a different proof. Note that at 
this point of the argument, <p may still depend on the sequence (t n ). We will prove below 
that this is not the case, and this implies convergence. We recall that the explicit formula 
for Pt(x,y) is 

Pf (x,y) = ^^e-I^^T/^ e -(^^)/^ sinh ^) . (3 ) 



It follows immediately that 



From the bound 
we conclude that 



P?(x,y) < Pt(x,y) 
2y d zd 



<p(.y, z) > 



(2n) d / 2 



In particular, for fixed z we find that cp(y, z) is a non negative harmonic function vanishing 
on the boundary and growing linearly with ya- This implies by the results of [B.] that the 
cone Vq, is of dimension two. 

Using the various symmetries and the convexity of the cone Vq we can say more about 
the structure of the function <p(x, y). First of all since Tn is of dimension two with u\ and 
U2 two extremals, and since <p(x,y) is harmonic in x at fixed y, non negative and zero on 
the boundary, we can find two functions fi and fa of y such that 



<p(x,y) = fi(y)u!(x) + h(y)u 2 (x) . 

Since u\ and u<i are not proportional, we conclude that fi(y) and f2(y) are harmonic, 
non negative and zero on the boundary of fi. Therefore writing fi and f 2 as nonnegative 
combinations of u\ and u 2 , and using also the symmetry of <p in x and y we obtain 

(p(x, y) = aui(x)ui(y) + bu 2 (x)u 2 (y) + c(ui(x)u 2 (y) + w 2 (x)wi(y)) , 

where a, b and c are three nonnegative constants which may depend on the sequence (t n ). 
Using the symmetry with respect to the hyperplane xa = 0, we conclude that a = b. 
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We now use the formula (see [C.M.SM.l]) for xaya > 

1 /■* 



= ^ X > V ) + j£faL (t-t)^ e-5/»<*-)e-l»-*|-/ a (*-.) p-(( f f o), y)d r . (5) 

If < the same identity without the term p^{x,y) holds. Notice that the function 
Ps((£;0),y) is symmetric with respect to the hyperplane yd = 0. Therefore if y* is the 
symmetric of y with respect to this hyperplane, we have for ya^d > 

Pt(x,y) =p^(x,y) +p t (x,y*) . 

Since u\(y*) = U2(y) and u 2 {y*) = u\{y) we conclude from (4) and our normalization of 
u\ (and U2) that a = 2(2n)~ d / 2 — c. It now remains to show that c = 0. 
We now assume xa > and split the integral 

x d f u I e-«3/2(*--) e -l*-€lV2(*-.) Pa(( ^o), y )df= J ff 1 + ^ 2j 



/o (t-«) 1+d/2 7D 

in two pieces given by 

Ki=x d £ /2 {t _ d s * 1+d/2 ^e-*2/ 2 («--)e-l^ a /2(*-.) pa((e -;o), y )df , 

^2 = * d T t + ^/2 / e-^^-^-i^T/ 2 (- S ) Ps(( ^o),,)^. 

./t/2 l J ~~ s ) 1 JD 

We now estimate each term separately starting with K 2 . Using Lemma 4 we have 
K 2 < r(l 4- «.(»)) £ tl+i/2 *_ t)1 ^ / o e-^««-)(l + ».(f.0))e-l*-« , « < -)<Je 

< r2 l +J /2 i -(l + J /2, (1+t , s(!/)) /■' / e -»3«'-»>(l +B ,(f,0)) e -l S! -«l 2 ««-»>df . 

it/2 (* ~~ S J ' </£> 

We integrate first over s and obtain 

^<0(l)rC^)(l + W .( I/ ))x, / ^±|4S> df. 

Jd K + IC-^I 2 ) 

Using the Herglotz representation for harmonic functions in a half hyperplane, we obtain 
finally 

K 2 < C(x)r^ +d ' 2 \l + v s (y))(l + o{x d )) , 
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where C(x) is bounded on compact sets in H d 1 . 
We now estimate K\ . We first observe that 

t^ 1+d ^x d f' 2 ds [ e-^ t e-\ s -^l t Ps {{lQ) 1 y)di<K l <{2n) d ' 2 Pt {x 1 y). 

JO J D 

Therefore, multiplying by t 1+d / 2 , taking t = t n and letting n tend to infinity, we obtain by 
the monotone convergence Theorem 

x d I ds p s {{£,0),y)d£ 
Jo Jd 

< (27v) d/2 (a(u 1 (x)u 1 (y) + u 2 (x)u 2 (y)) + c(ui(x)u 2 (y) + u 2 (x)u 1 (y))) . 
Dividing by Xd and then taking Xd to infinity, we obtain 

V(y) = ds p.((£o),y)d£< (2n) d / 2 (a + c)v s (y) . 

Jo J D 

In particular, this function is zero on the boundary of O. We now derive a better estimate 
on V. 

In equation (5) we take y = (fj, 0) and we integrate over ff and t. For any e > 0, we 
obtain i 

/ dt p t (x,(fj,0))dff 

J e JD , l^l^e" 1 

We now interchange the integrations in s and t and obtain 

/ dt p t (x,(ff,0))df] 
Je Jd , | ^7 1 <: e — 1 

= (SFpT 1 So ' ds (T^iy i + 3 7 I S D <if^' l ^ 2|2+l3 ' )/2( ' _s) f D ^ ?»_((£ o), (>?, o))dfi, 

^r*r^/^ e " (if "^ w / B ,«..^ £ - o) - ( * o)w 

- 0{l)Xd I * wi^r £ ' ds L pM ~ 0) ' ( " om ' 

Since V(y) is a non negative function (in fact harmonic) in the upper half hyperplane, we 
can let e tend to zero and using the monotone convergence Theorem we obtain 

V(x) < 0(l)x d [ dg V ^ 0) jj- < 0{l)x d [ dg ^ ( ^ 0) w/9 . 
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This is again the Herglotz representation and implies V(x) = C(x) (l + o(x d )), where C(x) 
is bounded on compact sets in R d_1 . 

We now come back to the estimation of K\. We have 

Ki < 2 i +d /2 t -(i+d/2) Xd f t/2 ds f ps (£ ),y)d£ 

JO J D 

< x d t^ 1+d ^O(l)V(y) < x d t^ 1+d ^C(y)(l + o{y d )) . 
We now take y d = —x d < large, x = y fixed. We have 

lim t 1 n +d/2 p tn (x,y) =(p(x,y) 

= 2au 1 (x)u 2 (x) + c(wi(» 2 + u 2 (x) 2 ) =cx 2 d + o(l + x 2 d ) < 0(l)x d (l + o{x d )) , 

which implies c=0. Therefore a = 2(2n)~ d / 2 for any accumulation point, which implies 
convergence and finishes the proof of the last part of Theorem 1. 

Finally, if t^ 1+d ^ 2 'pt(x, y) diverges for a pair of points (x, y) in by Harnack inequal- 
ity, it diverges for all pair of points. It follows from Lemma 4 that the cone Vn is of 
dimension one and the first part of Theorem 1 follows as in [C.M.SM.l]. This finishes the 
proof of Theorem 1. QED 

III. Asymptotic of the survival probability. 

In this section we give a proof of Theorem 2 on the asymptotic of the survival probability. 

We first consider the case where for some x £ O the function of time t 1 / 2 P x (T > t) has 
a finite accumulation point. In other words, there is a diverging sequence (t n ) such that 

tl/ 2 P x (T > t n ) has a finite limit. Applying Lemma 3, we conclude that tl^~ d ^ 2 p tn (x, x) 
is bounded. Therefore we are in the second case of Theorem 1, namely the cone Vq is of 
dimension 2. Applying Lemma 4 we conclude that t 1 / 2 P x (T > t) is bounded uniformly in 
t. On the other hand, any accumulation point is a non negative harmonic function which 
is symmetric with respect to the hyperplane x d = and zero on the boundary of O. This 
last statement follows (for t > 1) from Lemma 4 by the estimate 

VtP x (T>t)= [ Pl (x,y)ViP y (T>t-l)dy 
Jo, 



< Kyj j^-J^Piix, y) (1 + v s (y)) dy^K^j^-j (P X (T > 1) + v s (x)) . 

Therefore, any accumulation point is of the form Av s (x) where the constant A may depend 
on the accumulation point. We now show that this is not the case. Assuming x d > 
and integrating relation (5) with respect to y G O we obtain the following expansion for 
P*(T>t) 

P (T H > t) H Xd C — I e "( |2? "^ |2+ ^) /2(t " s) P - (T>s)d? 
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We also have easily from (3) 

We now split the integral in the above expression of P X (T > t) in two parts, namely 

Xd C — I e -\s-i\ 2 m-s) e -x 2 j2{t- S ) v _ (T s \ d f- Ml , Mo 

(27T) d / 2 J (t- S y+d/2 J D ^ [ >S)dt,-M 1+ M 2 

with 

Mi = — ' 
(2 

and 



Xd f I e -\3-t\ 2 /2(t-s) -xl/2{t-s) v _ ( T> s)d£ 



_ x d r* ds r w _$2 /2(t _ a) 2 /2{t _ s) _ 

~ (2n)*/* J t/2 (t-sy+d/2 J D e ^(€,0)^ > ■ 



We first estimate M\. We have, using Lemma 4, that M\ is bounded above by 

o (1)Xd r * . / e -Mvj«M / ig-gT+xgy /2 i + ,.«-o) - 

W Jo Mt-s)Jn V 2(< - s ' / (|2-H 2 + ^) / 

< o (1 ) Xd rf S , ( 1+V ' ll0) ,,A< OW Xd t-^ f 1 + V - il0) dl2 , 

- 7o V5(f - ») J D (|£ _ fi 2 + I -3) J/2 7 D (|^-«1 2 + *i) 7 

and since this last integral is proportional to the Herglotz integral, we obtain again 

M 1 <t- 1 / 2 C(x)(l + o(x d )), 
where C(x) is bounded on compact sets in R d_1 . For M 2 we have using again Lemma 4 

Integrating first over s, we obtain again the Herglotz integral and conclude as before that 

M 2 <0(l)t- 1 / 2 C(x)(l + o(x d )) . 
Finally we obtain for fixed x 

n f -l/2 

p x (t >t)= Xdl + 0((x d /Vif) + 0(l)C(x)r l / 2 {l + o(x d )) , 

V 27T 
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which implies A = y/2/ir. 

Assume now that t 1 / 2 P x (T > t) tends to infinity for some x G O. By Harnack 
inequality, this also holds for any y G O. Moreover from Lemma 4 we are in the case where 
the cone Vn is of dimension one. 

We now show that for any s6R and any x G O, we have 

P x (T>t + s) 

lim — — r — = 1 . 

t^oo P x {T>t) 

It is enough to consider s > since the case s < is obtained by taking the inverse of the 
above expression and replacing t by t + s. By the Gaussian bound and integrating first 
over y we have 



/ 

J V 



Pt{x,y)dy < I dy d = 



'\y d \<t 1/4 V*{T>t) 1/2 " ( 2nt ) 1/2 J\y d \<ti/*P x (T>ty/* "' V TV t 1 / 4 

This implies 

P (\X*\ < t^P x (T > if'* \T>t)< Jl — m , 

which tends to zero when t tends to infinity by our previous assumption. 
For s > the previous estimate implies 

w =E - (Pj - (T>i)|T>i) 

> E x (Px t (T > s) | |X t d | > t 1/4 P x (T > t) 1/2 , T > tj P x (|X t d | > t 1/4 P x (T > t) 1 / 2 | T > tj 

> E. (P X ,(T H > .) | |*<| > ( >/<p.(T > *)"* , T > ,) (l - Jl ^J^ 
From (3) we have 

P,(T H > s) = -== / U**-v*?/ 2 * _ e (- d + yd ) 2 /2 S \ dyd 

V27TS Jo ^ ' 



Therefore 



= 1 - -H= r e- 2 /2^ > x _ 2 v^ e -2/2- 
V27T Jid/v^ V27T a^d 

l> P a (T>t + g) > 



Px(T>t) 

_ -VtP x (T>t)/2s 



v 7 ^ (^P X (T > t)f' 2 )\ V tt (v ^ Px(T > £)) V2 

which tends to 1 when £ tends to infinity. 

For a fixed y G O, using reflection arguments, Harnack inequality and monotonicity 
in t to get away from the boundary, one can now show as in [C.M.SM.l] that as a function 
of x, any accumulation point of P X (T > t)/P y (T > t) is a non negative harmonic function 
which is zero at the boundary and equal to one for x = y. The conclusion of Theorem 2 
follows in that case since the cone is one dimensional. 
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Appendix: a symmetry estimate. 

Let n be a unit vector in R d_1 . Let y E with yd 7^ 0, we will assume for definiteness 
yd > 0, the argument is the same for yd < 0. Several quantities below will depend on y 
and n, but since we will not vary these two objects we will not mention this dependence 
explicitly. We define a domain 0+ by 

fi + = {x = (x, Xd) G O I \xd\ < (x — y) • n} . 

We also define two symmetry operators 5 ,± by 

S + (x) = (x + (xd — (x — y) • nj n, (x — y) • n) , 

and 

^"(x) = (f - (xd + (x - y) • n)n, —{x — y) • n) , 
Finally, we define the function q t (x,y) by 

=p t (S+ (»,?/) +pt(5~(a;),i/) -p t (x,y) . 
The following Lemma gives an estimate on p t (x,y). 
Lemma A. For any x G 0+ and any t > 0, we have qt(x, y) > 0, hence 

Pt(x,y) < p t (S + (x),y) +p t (S~(x),y) 

A similar result holds for the domain 

fi~ = {x E O | — \xd\ > (x — y) • n} . 

The proof is based on the maximum principle. We first observe that as a function of 
t and x, qt(x, y) is a solution of the heat equation. Moreover since yd > we have 

q (x, y) = 5(S + (x) -y) = S(x - S + (y)) . 

The boundary of Q + can be decomposed in three pieces. There is first the piece corre- 
sponding to Xd = (x — y) • n, where we have S + (x) = x, and therefore 

q t (x,y) =p t (S~(x),y) > . 

The second piece corresponds to Xd = —(x — y) • n, and on this hyperplane we have 
S~ (x) = x, therefore 

q t {x,y)=p t {S + {x),y) > . 

the last piece corresponds to 

D x {x d = 0} n {|x d | < (x-y) •fx] . 

On this set we have pt(x, y) = 0, and therefore 

q t (x,y) =pt(S + (x),y)+p t (S~(x),y) > . 

We now conclude the proof using the Phragmen-Lindelof version of the maximum principle 
(see [P.W.]). 
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